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Abstract. It has been suggested that a monitoring program employing the pixel lensing method to search for 
microlensing events towards galactic globular clusters may increase the statistics and discriminate among different 
halo models. Stimulated by this proposal, we evaluate an upper limit to the pixel lensing event rate for such a 
survey. Four different dark halo models have been considered changing both the flattening and the slope of the 
mass density profile. The lenses mass function has been modelled as a homogenous power -law for /j, £ {^i, ^u) and 
both the mass limits and the slope of the mass function have been varied to investigate their effect on the rate. 
The target globular clusters have been selected in order to minimize the disk contribution to the event rate. We 
find that a pixel lensing survey towards globular clusters is unable to discriminate among different halo models 
since the number of detectable events is too small to allow any reliable statistical analysis. 
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1. Introduction 

The detection of gravitational microlensing effect due to 
massive astrophysical compact halo objects (MACHOs) 
is undoubtedly one of the great success stories in 
astrophysics over the past decade. Since Paczynski's 
seminal paper (1986), microlensing surveys have begun to 
monitor millions of sources towards the Magellanic 
Clouds ULasserre et al. 2000" 'Alco ck et al. 2001|l . 
the galactic bulge (Udalski ct aL_1994^, the spiral 
arms (jPerue et ai. 200l')l and the Andromeda galaxy 
IjAuriere et al. 20011 ICalchi Novati et al. 20021 Pauhn- 
Henriksson et al. 2002a, b). The discovery of around 
20 candidates towards the LMC and more than a 
hundred events towards the bulge bears witness to the 
the potential of microlensing as a tool to investigate 
both the nature of halo dark matter and the structure 
of our Galaxy. Despite these successes a number of 
questions still remain. The optical depth measured 
towards the galactic bulge is at least a factor two 
larger than what can be accomodated by theoretical 
models (Bissantz ct al. 1997, Se venste r et al. 1999). 
The optical depth towards the LMC, measured by 
the MACHO group l| Alcock et al. 200 l|l . is too large 
by a factor ~ 5 to be accounted for by known stellar 
populations and too small by the same factor to be 
due to a halo full of MACHOs. Moreover, the implied 
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MACHO mass range, (0.1 — 1 Mq), is not easily 
reconciled with existing constraints on baryonic dark 
matter candidates (see, e.g., ICarr 199'4j) . even if it has 
been suggested that MACHOs could be not baryonic at 
all (Sazhin 1996, Gurcvich ct al. 19971). Furthermore, it 
is still not clear whether the lenses are indeed in the 
galactic halo since it is also possible that the contribution 
to the optical depth due to other, so far unknown 
stellar populations has been seriously underestimated 
HBennett 19981 IZhao 1999)) . In particular, the lensing 
rate due to stars or MACHOs within the LMC itself 
may be significant (|Auborg et al. 1999|ISalati et al. 1999L 
lEvans fc Kerins 20001 [Gyuck et al. 2000| ). 

However, it is important to stress that the inter- 
pretation of the results is clouded by our imperfect 
knowledge of the dark halo structure. In particular, 
studies show that it is not possible to constrain the 
halo structure and the lens mass function even with a 
high number of events if observations are taken in only 
one line of sight (Markovic fc Sommer - Larsen 1997)1 . To 
overcome this difficulty, the use of globular clusters has 
been suggested by many authors ( Gyuck fc Hol der 1997| 
IRhoads fc Malhotra 1998,Jetzer et al. 1998,) . This is pos- 
sible because globular clusters are distributed throughout 
the galactic halo, albeit concentrated towards the centre. 
Therefore, we can probe many lines of sight thus allowing 
a direct mapping of the variation of the density with 
position and a determination of the dark halo parame- 
ters such as its flattening. Furthermore, if the detailed 
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structure of the halo is lumpy, as suggested by some 
models of Galaxy formation, monitoring different lines 
of sight will average the results over such fluctuations 
in the halo density and thereby alleviate the potential 
worry about non- representative lines of sight. Finally, 
quick calculations ( |Gyuck fc Holder 1997| ) show that 
self - lensing may be for sure neglected in globular clusters 
so that the uncertainties about lens location are greatly 
reduced. Unfortunately, globular clusters are highly 
crowded (more than 100 stars arcsec^ at the core) thus 
leading to a small number of resolved stars. This means 
that classical microlensing (where one monitors resolved 
stars only) is not useful and pixel lensing IjCrotts 19921 
IBaillon et al. 1993|l must be considered. In this paper we 
estimate an upper limit to the event rate for pixel lensing 
observations towards globular clusters"'^ in order to inves- 
tigate whether a monitoring project could discriminate 
among different halo models. 

The plan of the paper is as follows. In Sect. 2 we briefly 
review some basics on pixel lensing in order to introduce 
the main formulae needed to estimate an upper limit to 
the predicted number of microlensing events. The model 
ingredients we need, regarding both the lenses and the 
targets, have been described in Sect. 3, while, in Sect. 4, we 
give our results for a selected sample of globular clusters. 
Conclusions are summarized in Sect. 5. 

2. Pixel lensing theoretical event rate 

Whilst in classical microlensing one monitors individual 
sources, in pixel lensing the sources are resolved only when 
they are lensed. We can therefore only monitor the flux 
in each detector element rather than the flux from indi- 
vidual sources. If a star is magnified sufficiently due to a 
lens passing close to the line of sight, then the total flux 
of the pixel containing the source star (coming from the 
lensed star, the other nearby unlensed sources and the sky 
background) will rise significantly above the noise level. 
The pixel lightcurve has then the classical Paczyiiski shape 
( ^Paczyhski 1986J as the usual microlensing lightcurve of a 
resolved star. But there is an important difference. Since 
the star emerges from the background only when it is 
lensed, the observed pixel lensing event timescale may not 
be a good tracer of the real Einstein time since we only 
see the portion of the lensed star lightcurve that emerges 
above the noise level. This means that the classical for- 
mulae used to evaluate the theoretical event rate and to 
estimate the optical depth from the data do not hold any- 
more. The first detailed study of pixel lensing was under- 
taken by Gould (1996) who defined two different regimes : 
a semiclassical regime, where the pixel flux is dominated 
by the source star and the observable timescale is a good 
tracer of the Einstein time, and a spike regime, where 

^ A similar study has been done by Rhoads & Malhotra 
(1998), but their approach was quite different from the one we 
will use here. We come back to this point in the conclusions. 



only high - magnification events can be detected, without 
any information on the Einstein timescale. Using Gould's 
formalism, Han (1996) predicted the pixel lensing event 
rate towards M31, but this has been obtained assuming 
a fixed sampling rate and unchanging observing condi- 
tions so that it is of limited applicability to ground - based 
observations. Here we do not employ Gould's formalism, 
but follow a different approach IjKerins et al. 200T)l which 
is more closely related to the pixel method. In order to 
better illustrate the approach we will use, it is helpful 
to briefly recall the basic steps of the analysis pipeline 
adopted for the lensing method (for a detailed description 
see lLe Du 2000IIGalchi Novati 2001|) . The starting point is 
a set of images of the same target taken on different nights 
and hence with photometric and seeing conditions chang- 
ing from one image to another. Before being corrected for 
seeing variations, the sequence of images must be geomet- 
rically and photometrically aligned to the same reference 
image so that the remaining variations may be due only to 
seeing effects or microlensing events and variable stars. To 
minimize the effect of seeing, one uses as base detector ele- 
ment an array of n pixels, which is termed superpixel. The 
number n is fixed in such a way that the superpixel an- 
gular dimensions are larger than the median seeing of the 
images. For instance, if the median seeing is l.Sarcsec and 
the pixel size is O.SOarcsec, we choose n = 6 so that the 
superpixel size is l.Sarcsec, larger than the median see- 
ing. Binning the fiux in the superpixel reduces the effect 
of seeing variations, but this is not enough to make them 
negligible. To this aim, one turns to the pixel method, an 
empirically derived statistical correction which is applied 
to each image to match it to the observing characteristics 
of the reference image^. 

Let us consider now a superpixel collecting the pho- 
tons emitted from a set of stars in the superpixel field of 
view. Suppose that an unresolved source is undergoing a 
microlensing event and let us call the source magnifica- 
tion factor A{t). After alignment and seeing corrections, 
the excess superpixel photon count AA'gpx on an image i 
obtained at epoch ti due to a microlensing event is : 

AN.p^iU) = Nki[A,p^{U) - 1] = NsU,[A{U) - 1] , (1) 

where Ng and A'bi are the source and baseline photon 
counts in the absence of lensing and /soo the fraction of 

^ The pixel method corrects for the different seeing of the im- 
ages by linearly rescaling the flux of a given superpixel on each 
image. The scaling coefficients are determined (on an image - 
by - image basis) through a statistical analysis of the differences 
among the median of the current image and that of the refer- 
ence image. As a result, the superpixel lightcurve turns out to 
be ffat (within the uncertainties) unless a microlensed star or 
a variable star is in the field of view of that superpixel. The 
pixel method has been successfully implemented and has been 
proven to be an efficient tool to look for microlensing events to- 
wards the Andromeda galaxy (Ansari et al. 1997, 1999, Auriere 
et al. 2001, Calchi Novati et al. 2001, Paulin - Henriksson et al. 
2002a,b). For further details, we refer the reader to the litera- 
ture (see, e.g., [Kaplan 19981 ITelH 2000.. Calchi Novati 20'on ) 
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the seeing disc contained within the superpixeP. It is im- 
portant to stress that we are defining ^spx as the super- 
pixel count variation factor which acts as the observable 
counterpart of A since the only quantities which may be 
detected from observations are /see-^s(^ — 1) and A'^bi, but 
not Ns and A separately. If we denote the unlensed appar- 
ent source magnitude by nis and the surface brightness (in 
mag/arcsec~^) of the target galaxy and of the sky respec- 
tively as /Zgai(a;, y) and /isky, and iVbi can be estimated 
as : 



A^bi = No rispxx 10 — + 10'— 



(2) 
(3) 



where rJgpx is the superpixel angular area and {x, y) its 
coordinates in a rectangular system with origin on the tar- 
get centre and oriented along its main axes. Note that we 
are assuming that both /igai and /igky are constant along 
the full sequence of images. Because of changing observing 
conditions, this could not be true, but we may assume that 
the photometrical aligment has reduced the fluctuations 
in /isky to a negligible level. 

The superpixel noise on image i is : 



or, Si^J Nspx{ti) 



where : 



(5) 



The threshold noise level ax is determined by the super- 
pixel flux stability, while the correction factor 6i takes ac- 
count of the fact that the pixel method is not photon noise 
limited. A signal is regarded as being statistically signifi- 
cant if : 

ANsp^iU) > 3(7, . 

This means that a microlensed source is detected if the 
magnification is higher than a minimum value which may 
be computed from Eq. as : 



^min(^j) — 1 + 



3a,, 



fscc-Ns 



(6) 



A special case occurs when ai = ut, giving a threshold 
magnification l|Kerins et al. 2flOT| : 



At = 1 + 0.0075 



bl 



/sec N g 



r^snx 10 2.5 -(- 10 2.5 

1 + 0.0075 

/sec 



(7) 



10" — 

^ This may be easily estimated as /sec ~ ""'"psf / ^^spx with rpsf 
the radius of the seeing disk (in arcsec). Note that /sec < 1 
given the way we choose the superpixel dimensions. 



Eq.© may be used to define a threshold dimensionless 
impact parameter ut- The relation between A and u is 
the same as in classical microlensing ( |Paczynski 1986| l : 



A 



(8) 



Inverting this equation to get u — u{A) and using Eq.Q 
we can obtain the value of the impact threshold parameter 
as : 



Ut 



V2 



An 



1/2 



(9) 



The dependence of ut on Ng means that the pixel event 
rate depends on the source luminosity function 4>{M) with 
M the absolute magnitude, while the dependence on A^bi 
leads to a spatial variation of ut- Note that ut depends on 
the observing filter too. Although in Eq.@ Nq is dropped 
out, the dependence on the filter is implicit in the luminos- 
ity function since this latter quantity changes according to 
the observing photometric band. 

We can now compute a theoretical pixel lensing event 
rate at sky coordinates (x, y) as : 



Tpi = {uT){x,y) rc{x,y) , 



(4) where we have defined : 



{uT){x,y) 



J UT{M,x,y) <j){M)dM 
/ (t>{M)dM 



(10) 



(11) 



In Eg. ljlOII . Tcix^y) is the classical microlensing event 
rate which we can estimate as (|R.oulet fc Mollerach 19971 
l^tzer 1998a : 



^c{x,y) = 



2 T{x,y) 



TT {t^){x,y) 
where the optical depth r is given by : 

p[s) s[l — s)as 



(12) 



(13) 



where Dg is the distance of the target (assumed to be 
within the galactic dark halo) and p the halo mass den- 
sity expressed as function of s = Dd/Dg (with the 
distance to the lens). In Ea. (|12|l . {ts) is the Einstein ra- 
dius crossing time averaged over the lens mass function 
and the halo transverse velocity distribution. We use the 
same procedure described in Han & Gould (1996) to es- 
timate this latter quantity, but we do not report here the 
explicit expressions for the sake of shortness. 

The pixel lensing theoretical event rate Fpi estimated 
by Ea. (|10(l must be considered as a strict upper limit to 
the really observed event rate. Actually, Ea. HlOl) assumes 
that the observing conditions are unchanged which is not 
exactly true even if the pixel method partially reduces the 
variations due to seeing effects and atmospheric condi- 
tions. Besides, since one usually has ^min > At, Eq.l^UJ 
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also tends to overestimate the true mean pixel lensing 
cross section. Finally, Ea. ((T^ assumes perfect detection 
efficiency, whilst the efficiency function for pixel lensing 
surveys is not known at all (but surely much lower than 
one). These considerations suggest that Fpi cannot be di- 
rectly compared to observations. Nonetheless, Ec ljlUI) is 
enough for our preliminary study. Actually, we want to 
estimate the predicted number of events towards globu- 
lar cluster to see whether a pixel lensing survey towards 
these targets may be worthwhile or not and whether the 
information one can gain is useful in discriminating among 
halo models. 

3. Model ingredients 

There are different physical ingredients entering Eq.lQUJ 
that we have to fix in order to estimate the theoretical 
pixel lensing event rate towards globular clusters. One 
needs to assume a dark halo mass density profile to com- 
pute the optical depth r with Eq. (|13|l , while also the trans- 
verse velocity distribution and the lens mass function enter 
the estimate of (i^). On the other hand, the source prop- 
erties, namely the surface brightness and luminosity func- 
tion, have to be assigned to estimate the average thresh- 
old impact parameter (ut). In this section, we describe 
our choices for each of these ingredients. As a prelimi- 
nary remark, we stress that our main aim is to investigate 
whether a survey looking at globular clusters may be used 
to discriminate among dark halo models. To this aim, it is 
useful to observe globular clusters which are far above the 
galactic plane in order to minimize the disk contribution 
to the lensing rate. The criteria we employ to select our 
main targets will be described in Sect. 4. For the time be- 
ing, let us stress that the only galactic component where 
we assume the lenses could be is the dark halo so that we 
have to model only this galactic component. 

3.1. The dark halo 

Modeling the dark halo is not an easy task given the 
paucity of information we have on the dark matter dis- 
tribution. Hence, there are different models which fit 
equally well the same constraints. Furthermore, our po- 
sition within the Galaxy complicates the measurement of 
rotation curve which usually gives the strongest constraint 
on the dark halo, when studying external galaxies. We ex- 
plore here four different halo models changing both the 
shape and the slope of the mass density. 

The first class of models we consider is the non singular 
isothermal spheroid with mass density : 



p(i?,z) 



Pq 



Rl 



Rr- 



q arcsin (-^/l — q) R"^ + / q^ + R^ 



(14) 



where is the local dark matter density, q the halo flat- 
tening, i?o the distance of the Sun to the galactic centre 
and Rc the core radius. The local parameters are fixed as : 

Pq = 0.008 Mq/pc'^ , i?o = 8.0 kpc , 



whilst we explore models with three different values of 
q = 0.5, 0.75, 1.0. Finally, the halo core radius is fixed by 
the following relation : 



I" ( oo ) arcsin ( vT~~g) 



Rc 



i?5 



1/2 



(15) 



where Uc(oo) is the asymptotic value of the flat rotation 
curve flxed to be 200km/s. The mass density is all we 
need to evaluate the optical depth towards a target with 
galactocentric coordinates {Dg, l,b) : 

I h) 47rGJ^ jR'o + Rl) ^ ^ 



q arcsin (-^/l — q) 



s(l-s) 



A,. 



BisoS + Cii 



ds 



(16) 



with : 



Aso = Dl (cos^ b - 



-2 • 2 

q sm 



Biso — — 2i?o-Ds cos 6 cos Z , 

Ciso — Rq + -Rc ■ 

To estimate the average time duration (ts), we need to 
know the transverse velocity distribution of the lenses. 
This should be obtained in a dynamically self- consistent 
way starting from the mass density to recover the model 
distribution function in the conflguration space and then 
integrating over the spatial coordinates and the velocity 
parallel to the line of sight. This is not a difficult task for 
spherical models (see, e.g., |Binney &: Tremaine 1987| l, but 
the result could be complicated to obtain for spheroidal 
models. However, we can follow the usual approach as- 
suming that the transverse velocity distribution for the 
isothermal models is isotropic and maxwellian : 



2vi 



-'H 



exp 



dv± 



(17) 



with vh = 210 km/s. 

The second model we consider is the so called 
truncation flat model (hereinafter TF model) which 
has been successfully used to reproduce the dynam- 
ics of the Magellanic Stream IjLin et al. 1995|l and 
the motions of the Milky Way's satellite galaxies 
(Wi lkinson & Evans 1999)1 . The mass density is spheri- 
cally symmetric and it is given as : 



p{r) 



(18) 



where Mtot the total mass of the halo and a scale radius. 
We fix these parameters as found by Wilkinson & Evans 
(1999): 



tot 



1.9x10^2 M, 







170 kpc 
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Inserting Eq. (|18|l into Eq. H13|) we get the following expres- 
sion for the optical depth : 



s(l-s) 



P{s) [P(s) + r2]3/2 



ds (19) 



with : 

P{s) = Atfs"^ + Btfs + Ctf 
with : 



Atf = D^^ , Btf = -2RqD, cos & cos ; , Ctf = Rl ■ 

Wilkinson & Evans (1999) have computed the distribu- 
tion function for the TF model. We have numerically in- 
tegrated their expression in the case of isotropy in the 
velocity space and interpolated the results to get the fol- 
lowing approximated expression for the transverse velocity 
distribution of the model we are using : 



f{v±) dv± cx vj_ exp 



with : 



4.151 



0.5984 Vs 



1.0379 



dv^{20) 



V, = y/GMtot/rs 219 km/s . 

The last ingredient we need to fully characterize our halo 
models is the local mass function (hereinafter MF) of the 
lenses*. In the usual approach, one assumes that all the 
MACHOs have the same mass so that the MF is a 5 - Dirac 
centered on the common mass /i (expressed in solar units). 
However, this distribution is not realistic and can lead to 
seriously wrong results since the average time duration is 
quite sensitive to the MF. A first step towards a more 
realistic approach is to use a power - law MF : 



(21) 



The slope a and the lower and upper mass limits of the 
lenses (/i;,/x„) are poorly constrained. However, compar- 
ing the data towards LMC with the theoretical expec- 
tations (also taking into account the detection efficiency 
and the self- lensing effects), Cardone et al. (2001) have 
estimated the slope a of the local MF using the same 
spheroidal models we are employing here. Even if the 
model parameters used are slightly different, we will use 
their results in our computations. They did not consider 
the TF model in their analysis so that we will fix the 
parameters of the MF for the TF model to the values 
they have obtained for their isothermal spherical mod- 
els. Although quite arbitrary, this choice is enough for 
our aims and more realistic than the usual common mass 
hypothesis. Furthermore, we will consider two values for 
/i; = 0.001, 0.01, whilst the upper mass limit is fixed to be 
flu = 1-0 IjCardone et al. 200111 . Table 1 summarizes the 
main properties of the different models considered and in- 
troduce the labelling we will use in the following. 

* Following Han & Gould (1996), we are implicitely assuming 
that the MF is homogenous, i.e. dn/dfi = p{R,z)/ pQX^{fi), 
so that we have only to give the local MF ip{p) (see also 
IGrenacher et al. 19991 lUardone et al. 2001t . 



Id 


p{R,z) 
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QO 


(«!, Qii) 


I50a 


Eq.ini 


Eq.ini 


0.50 


1.10 


(0.50, 1.48) 


I75a 


Eq.ini 


Eg.lTTl 


0.75 


1.17 


(0.80, 1.50) 


IlOOa 


Eg.lini 


Eg.llTTl 


1.00 


1.21 


(0.84, 1.56) 


TFa 


Eq.lUHl 


Eg.CT 


1.00 


1.21 


(0.84, 1.56) 


I50b 


Eq.linj 


Eg.jlVl 


0.50 


1.12 


(0.27, 1.69) 


I75b 


Ea.dTIl 


Eq.ini 


0.75 


1.34 


(0.77, 1.89) 


IlOOb 


Eq.CU 


Eq.im 


1.00 


1.34 


(0.77, 1.90) 


TFb 


Eq.lUHl 


Eg.lCT 


1.00 


1.34 


(0.77, 1.90) 



Table 1. Labelling and main properties of the models 
considered in our analysis. The first half of the table refers 
to models with {fii,fiu) = (0.001,1.0), whilst the second 
one is for models with (iii^Hu) = (0.01,1.0). The lower 
and upper limit of the range used for a are (a;, Uu) with 
ao as best fit value l|Cardone et al. 20111)1 . 



3.2. The globular cluster model 

In order to use Eq. l(TT))l to evaluate the pixel lensing theo- 
retical event rate, we need to give also the surface bright- 
ness of the target. The globular cluster surface brightness 
is very well fitted by the King model (King 1962, 1966) : 



mc{R) = A*o + 5 log 



5 log 



1 - 



(l + r2/r2)V2 



(l+i?2/r2)l/2 (1+^2/^2)1/2 



(22) 



with Ho the central surface brightness (in mag/arcsec^), 
rt the tidal radius and the core radius. These 
parameters will be fixed as reported in the cat- 
alogue^ of galactic globular clusters compiled 
by Harris (1996) which can be found under 
|http : //physun . physics . mcmaster . ca/Globuleir . htinl| 
A key ingredient in the evaluation of the event rate by 
Eq. (|1U|I is the target luminosity function (hereinafter LF) 
which enters the determination of the average threshold 
impact parameter {ut). However, the measurement of 
the globular cluster LF is a quite complicated task 
because of the severe crowding. Typically, the use of 
the Hubble Space Telescope allows one to resolve a 
sufficient number of stars, but the brightest ones cannot 
be observed because of saturation of the WFPC instru- 
ment. A possible solution is to combine ground -based 
observations of the brightest (and resolved) stars in the 
cluster with the faintest star observed with the HST. 
Unfortunately, this method has been used to determine 
the LF only for four of the globular clusters we have se- 
lected, namely NGC52 72 IIKood et al. 199911 . NGC6254 
fPiotto fc Zoccah 1999|l . NGC6341 and NGC 7078 
(Piotto et al. 1997J. Furthermore, the measured LF is 
incomplete at the lower end, while a correct estimation 
of {ut) invokes the knowledge of the LF over the full 
brightness range. Thus, we choose to adopt a LF derived 



We have used the 1999 version. 
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via numerical simulations rather than an experimental 
one. The use of theoretically derived LF allows us to 
take into account the presence of stars in all evolutionary 
stages: from the faint and richly populated end of the 
Main Sequence, to the bright and fast Asymptotic Giant 
Branch (AGB) and Post-AGB phases. The LF used 
in this work have been calculated implementing the 
evolutionary code developed for the synthesis of stellar 
population by Brocato et al. (1999, 2000). The results 
of this code have been largely tested on the young LMC 
clusters, and on the old galactic globular clusters showing 
a very good agreement. 

We briefly report here the main assumptions and in- 
gredients of the code. All the computations refer to simple 
stellar populations, i.e. single -burst stellar systems with 
homogeneous chemical composition. Once the age T and 
the chemical composition {Z, Y) of the stellar population 
are fixed, the mass spectrum of the evolving stars in the 
population is obtained assuming a Scalo - like Initial Mass 
Function, randomly populated through a Monte Carlo 
method. The physical quantities associated to each evolv- 
ing mass (luminosity, effective temperature, surface grav- 
ity, etc.) are derived adopting the scenario of stellar evo- 
lution of the Teramo- Pisa- Roma group (see Cassisi et 
al. 1998 and references therein). Hence, to relate the the- 
oretical quantities logL/L© and logTefj with the obser- 
vational colours and magnitudes, we assumed the atmo- 
spheres models by Lejeune et al. (1997). For a detailed 
description of the complete computational procedures of 
the code we remind to the quoted Brocato et al. papers. 

Being interested to the galactic globular clusters, in the 
following we refer to LF obtained from simulations of sim- 
ple stellar populations with a fiducial age T = 15 Gyr and 
for three different chemical compositions: Z = 10^"', Y ~ 
0.23 (metal- poor populations); Z = 10^^,1" = 0.23; 
Z = 10^^, Y = 0.25 (metal-rich populations). 

4. Selection of globular clusters and results 

The catalogue of globular clusters compiled by Harris 
(1996) contains almost 150 systems and many of them are 
distributed in the galactic halo. In order to best probe the 
dark halo structure, it is important to reduce as much as 
possible the contribution of the other galactic components 
to the event rate so that we can safely assume that the 
events that will be (eventually) detected are due to halo 
lenses. To this aim, we will consider as targets only those 
globular clusters which are far above the disk plane, i.e. 
with |6| > 10°. This constraint reduces also the bulge con- 
tribution. However, to completely exclude the possibility 
that the line of sight goes through the bulge, we further 
impose X < 6kpc, where X is the a; coordinate of the 
cluster in a rectangular system centered on the Sun with 
the X - axis positively directed towards the galactic cen- 
tre. We also ask that our targets reside inside the galactic 
halo so that we will not consider globular clusters with 
Dqc > 100 kpc, with Dqc the distance of the cluster to 
the galactic centre. Using these selection criteria we have 



Name 


TI30 


T/75 


T/lOO 


ttf 


NGC 288 


1.22 


1.06 


0.91 


1.48 


NGC 1904 


2.27 


1.69 


1.34 


1.95 


NGC 2419 


12.3 


9.74 


7.86 


9.48 


NGC 4147 


3.30 


3.39 


3.31 


4.96 


NGC 5024 


3.23 


3.38 


3.36 


5.17 


NGC 5053 


2.93 


3.01 


2.95 


4.61 


NGC 5272 


1.75 


1.61 


1.46 


2.38 


NGC 5466 


3.14 


3.23 


3.18 


4.99 


NGC 5904 


1.61 


1.39 


1.23 


2.56 


NGC 6171 


1.66 


1.31 


1.09 


2.97 


NGC 6205 


1.51 


1.22 


1.02 


1.83 


NGC 6218 


0.89 


0.67 


0.55 


1.22 


NGC 6254 


0.72 


0.54 


0.44 


0.95 


NGC 6341 


1.67 


1.31 


1.09 


1.87 


NGC 6366 


0.48 


0.35 


0.28 


0.58 


NGC 7078 


2.64 


2.07 


1.70 


2.97 


NGC 7089 


3.06 


2.63 


2.27 


4.10 


Pal 13 


5.79 


5.38 


4.81 


6.94 



Table 2. Optical depth (in units of 10~^) towards the 
globular cluster centre for the selected targets. Note that 
we have dropped the label a or b from the model id since 
T does not depend on the MF. 



compiled a sample of potentially interesting targets among 
which we only select those clusters which are observable 
from the TTl telescope^. For these systems, we have eval- 
uated the theoretical pixel lensing event rate using Eq. (|10ll 
assuming that MACHOs made up 20% of the halo dark 
matter l|Lasserre et al. 200(11 lAlcock et al. 200i| . 

Table 2 reports the optical depth towards the centre 
of the selected globular clusters and shows that r is best 
suited to discriminate among different halo models since 
the values predicted are quite sensitive to both the flat- 
tening and the slope of the density profile. Furthermore, r 
does not depend on the MF and hence possible degenera- 
cies between the density parameters and the MF slope are 
avoided. Unfortunately, measuring the optical depth from 
pixel lensing events is not immediately possible. The usual 
estimator, r,„cas = {t^ l'2'E)^tE,i/ s{tE,i),'ai&y'i^ot be used 
since the Einstein time is is not available from pixel lens- 
ing events. An alternative estimator has been proposed 
IIGondolo 1999|l . but it turns out to depend on the lumi- 
nosity function and also requires a detailed knowledge of 
the efficiency of the pixel lensing survey. Unfortunately, 
both these quantities are difficult to obtain or not avail- 
able at all. 



The TTl telescope is situated in Toppo di Castelgrande 
(Potenza, Italy) with coordinates (40.6° N, 15.8° E) 
and will be used by the SLOTT- AGAPE collabo- 
ration to monitor the Andromeda galaxy searching 
for microlensing events by the pixel lensing method. 
Details on the TTl may be found at the web site 
http: //www. na. astro . it/oacmedia/ttl_site/index .htmll 
whilst further informations on the SLOTT -AGAPE project 
are available at , http: //www. na. astro . it/ slag/ slott .html, 
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Name 


rc(/50a) 


rc(/75a) 


rc(/100a) 


r^TFa) 


rc(/50b) 


rc(/756) 


rc(/ioofe) 


rc(TF&) 


NGC 288 


2.4li-t 


2 1+^-^ 






1 6+"''* 




1-31^1 


2.4li-^ 


NGC 1904 


3.6+?-? 


2.8+i-ij 


2.4t?-; 


4.1+?1 


2.5 Q g 


2.0+i-? 


i-elo'o 


2.7+1-^ 


NGC 2419 


10.1+i'^ 




7.2+^-3 


10.8+^'t 


6.8+1^ 


5.9+1-'^ 


4-7+?-8 


7.0+11 


NGC 4147 


4.6l|i 




4.9+^-'3 




3.111-1 


3.3111 


3.2111 


5.6+1? 


NGC 5024 


4.6111 




5 1+4:5 




s.itii 


3.4111 


3.4111 


6.0111 


NGC 5053 


4.0l?1 




4.5+3-? 




2.7tli 


3.0+1-f 


2.9+1-f 




NGC 5272 


s.otli 


2.9+\i 


2-8+?1 






2.0lJi 


1.81J-? 


3.511-8 


NGC 5466 


4.51^-^ 




4.9l^-^3 




3.1111 


3.3111 


3.2111 


5.9111 


NGC 5904 




2.9+li 


2-8+?1 


6.6+1-1 


9 9+1.2 
—0 8 


2 1+1-^ 

— U 8 


l-8l°o1 


4.412-3 


NGC 6171 


3.6+17 


2.9+li 


2.7+f^ 


8.6tli 




2.llcV8 


1 7+0:9 

-L • ' _ (5 5 


5.61^1 


NGC 6205 




2.5tl't 








1-81E1;; 


i-5i;if, 


3.111-? 


NGC 6218 


2.2+i-? 


1 7+0.9 

-L • ' _ (] ^ 






i-5ir6 


1.2l«-6 


i.oi;i'^ 


2.611-1 


NGC 6254 


1 9+^-1 

— 9 


1 C+O.'S 


1 0+1:1 

— 6 




i.3i;;l 


i-ol?,1 


o.8i;it 


2 1 + ^-1 

— 8 


NGC 6341 


o 0+2.0 
"J"J-1.6 


2 7+^-* 

^- ' -1.2 


2 4+2'i 




9 9+1.2 
^•^-0.8 


-"-••^-0.7 


-'-•O-0.6 


1+1:6 


NGC 6366 


1.41°:? 






2.2i-;;f, 


o.giEJi 


7+"-" 
' -0.3 


U.D_o 2 


1 4+0-7 

-^•^-0.5 


NGC 7078 




3.7+?-? 




6 7+''''^ 
' -3.1 


3.111;^ 


«+1.3 
^•"-1.0 


2 1 +-'-^ 
^•-^-0.8 


4.4l?j 


NGC 7089 




4.51^;;; 


^•-■--2.0 


7+7.6 

^'•'-4.0 


3.411:^3 


9+1.6 
o.^_i 2 


2.711* 


t; 7+3.0 


Pal 13 


7 1+4.3 
' •-'--3.4 




6 2+^-'' 
"•^-2.9 


10.7+^-^ 


4.8ll:« 


4.61^7 


4.ll?l, 


7 0+'^-^ 

' •'-'-2.6 


Total 




62 2+-'^-^ 


61 1+-''^ 


119.01^^3" 


48.9lf3-.g 


45.6l-;^l 


40.11™- 


70 1 +40.7 
'8.1_29.0 



Table 3. Theoretical classical microlensing event rate (in units of 10~^ events/star/year) towards the centre of the 
selected globular clusters. For each model, we report the quantity XXlJJ meaning that the estimated range for Tc 
is {XX — ZZ, XX + YY) obtained by varying the slope a of the lensed MF in the interval (a;, a„), -while XX is the 
value of Fc for the MF -with slope qq- Both (a;, a„) and ao are given in Table 1. 



Tabic 3 gives the classical event rate to-wards the cen- 
tre of the selected globular clusters. The dependence of Fc 
both on the halo models and the MF is clearly evident. 
Actually, an accurate estimate of Fc is not possible unless 
the uncertainty on the slope a of the lenses MF is signi- 
ficatively reduced. Ho-wever, even if a were known with 
much less uncertainty, observations towards a single glob- 
ular cluster can not uniquely discriminate among differ- 
ent halo models. For instance, Table 3 shows that the two 
extreme models I50a and TFb predicts almost the same 
value (within the uncertainties) of Fc towards all clusters. 
It is, however, interesting to observe that the quantity 
Tl/rl"*, with ri the rate towards the i-th cluster and 
r*°* = X^jF*, is almost independent on the MF and can 
be used to discriminate among different halo models. The 
same results approximately applies if we use Fpi instead 
of Fc. Since Fpi can be easily measured from observations, 
this could be an interesting method to discriminate among 
different halo models avoiding the uncertainties related to 
the lenses MF and the need for the measurement of r. 

While T and Fc do not depend on the target LF, this 
latter quantity plays a key role in the determination of 
{uT)(x,y). In Table 4, we report the main photometric 
parameters for the sample of globular clusters we are con- 
sidering and the value of {ut) defined as follows : 

/ {uT){x,y) X io-A'Gc(x,2,)/2.5 dy 



Name 


Mo,n,rc 


Hv 


Mv 


Z 


{ut) 


NGC 288 


14.43, 12.94, 1.42 


14 


69 


-6 


60 


10 





27 


NGC 1904 


16.23, 8.34, 0.16 


15 


59 


-7 


86 


5 





26 


NGC 2419 


19.83, 8.74, 0.35 


19 


97 


-9 


58 


1 





01 


NGC 4147 


17.63, 6.31, 0.10 


16 


48 


-6 


16 


3 





13 


NGC 5024 


17.39, 


21.75, 


0.36 


16 


38 


-8 


77 


2 





14 


NGC 5053 


22.19, 


13.57, 


1.98 


16 


19 


-6 


72 


1 





10 


NGC 5272 


16.34, 


38.19, 


0.55 


15 


12 


-8 


93 


5 





22 


NGC 5466 


21.28, 


34.24, 


1.64 


16 


15 


-7 


11 


1 





10 


NGC 5904 


20.49, 


12.07, 


1.96 


14 


46 


-8 


81 


10 





33 


NGC 6171 


18.84, 


17.44, 


0.54 


15 


06 


-7 


13 


20 





19 


NGC 6205 


16.80, 


25.18, 


0.78 


14 


48 


-8 


70 


5 





26 


NGC 6218 


18.17, 


17.60, 


0.72 


14 


02 


-7 


32 


6 





34 


NGC 6254 


17.69, 


21.48, 


0.86 


14 


08 


-7 


48 


6 





33 


NGC 6341 


15.58, 


15.17, 


0.23 


14 


64 


-8 


20 


1 





29 


NGC 6366 


21.24, 


15.20, 


1.83 


14 


97 


-5 


77 


30 





26 


NGC 7078 


14.21, 


21.50, 


0.07 


15 


37 


-9 


17 


1 





27 


NGC 7089 


15.92, 


21.45, 


0.34 


15 


49 


-9 


02 


4 





22 


Pal 13 


23.36 


2.20, 0.48 


17 


31 


-3 


51 


4 





05 



Table 4. Main properties of the sample of globular clus- 
ters selected as described in the text. The King model 
parameters are the central surface brightness /iq (in the 
V- mag/arcsec^), the tidal radius rj and the core radius 
Tc (in arcsec). The distance modulus and the inte- 
grated absolute visual magnitude My fix the luminosity 
of the cluster. The last two columns gives the metallicity 
Z (in units of 10~^) and the value of (ut). 



{ut) 



J io-/'Gc(x,y)/2.5 dx dy 



(23) 



which is the average threshold impact parameter spatially 
weighted with the cluster surface density. We will use this 



quantity later to approximatively estimate the number of 
observable lensing events. 
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5. Conclusions 

The results in the previous section have shown that ob- 
servations towards globular clusters have the potential to 
discriminate among halo models and to shed light on the 
MF of the lenses. Actually, in order to see whether a pixel 
lensing survey may indeed be useful, we must evaluate 
(at least qualitatively) the theoretical number of events 
predicted for a given model. We can estimate it as : 

7Vp7 ~ iV4m < mii„0 X T„nobs x re(0, 0) x (24) 

where Ni,{m < miim) is the number of cluster stars 
brighter than and Tobs the total period of obser- 

vations, that is the product of the observability^ per year 
and the duration of the monitoring campaign. In Eg . (1241) . 
we have approximated : 



rc(x,?;) ~ rc(0,0) , {uT){x,y) ~ {ut) ■ 

While the first approximation works very well (within 
2%), the second one leads to a slight overestimation of 
the number of events, but nonetheless is useful to get or- 
der of magnitudes results that are enough for our aims. 

To estimate Af*(m < mii,„), we have to choose a value 
for the limiting magnitude. This is a complicated task. 
In principle, all the stars in the cluster could be sources 
of a detectable pixel lensing event provided that the am- 
plification is sufficiently high. Actually, it is highly un- 
likely that very faint stars give rise to detectable events. 
Furthermore, the limiting magnitude depends on the effi- 
ciency of the monitoring campaign and of the data anal- 
ysis pipeline. We fix niu^ — 26 for NGC 7078 and scale 
miiiii to the other globular clusters according to the dis- 
tance modulus. A higher value of rniim should reduce the 
predicted number of events. To estimate 7V^(to < mii,„), 
we have first to renormalize each simulated LF using the 
condition : 

Cnorm J lO' (j){M)dM = (25) 

with My the integrated absolute magnitude of the cluster. 
Then, the number of stars brighter than miim is simply : 

7V,(m < mii„) = Cnorm / ^{M)dM (26) 

where Cnorm is the normalization coefficient, Afmin the 
minimum value of the absolute magnitude and /iy the 
distance modulus. The total number of events predicted 
for pixel lensing observations is reported in Table 5 for 
a monitoring campaign lasting three years. To study the 
impact of the metallicity of the cluster population on 
the estimated number of events, we arbitrarily assume 

The observability is strongly dependent on the telescope 
site. We have assumed that observations are performed using 
the TTl telescope and deem a cluster as observable if it is 
higher than 20° on the horizon for more than 3 hours per night. 



Model 


Z = 10"^ 


Z = 10^^ 


Z = lO^"* 


I50a 
I75a 
IlOOa 
TFa 
I50b 
I75b 
IlOOb 
ITFb 


''J-0.40 
n 70+0-55 
"-'-0.32 

0.66t°:^3l 

^••^"-'-0.63 

0.49«:?^ 
0.46«:?^ 
0.43«;?^ 


1.00_Q 74 

1 29+^°° 
1 20+^"'^ 

^•^'-'-0.56 

o.9ot°:i 

0.86«:S 

' ''-0.29 

1 63+°**'' 

^•"■J-0.61 


1 70+ -^■'''^ 

1 f.n + 1.26 
1 1-9+1.33 

14+2.75 

1 -iq+0.68 
J^-^'J-0.48 

^^■^"-0.38 

^^•""-0.37 

2 06+^ °^ 



Table 5. Total number of events for the different models 
considered and the three simulated luminosity functions 
(with Z the metallicity of the population). The central 
value refers to the MF with slope ao, while the lower and 
upper limits have been obtained varying the slope a in the 
range given for each model in Table 1 . 



that all the clusters of the sample have the same metal- 
licity and use simulated luminosity functions for stellar 
populations having three different values of Z, as de- 
scribed before in Sect. 3. The observed trend with Z is 
easy to explain qualitatively. The lower the metallicity, 
the fainter the stars in the same evolutionary phase and 
hence the higher the threshold amplification At needed. 
Since oc ut oc A^^, it is expected that N^^ decreases 
as the metallicity Z increases. 

Unfortunately, whatever is the value of Z, it turns out 
that the number of expected events is so small that all the 
tests we have suggested to discriminate among halo mod- 
els are useless because almost no events are predicted. It 
is important to remember that the results in Table 5 have 
been obtained using Ea. HlUI) which assumes a perfect de- 
tection efficiency and ideal observational conditions. Thus, 
the numbers in Table 5 should be considered as strict up- 
per limits, the actual number of detectable events being 
far lower. 

Note also that, up to now, all the quantities reported 
have been given without any uncertainties. Actually, both 
the halo model parameters and the globular clusters prop- 
erties (King model parameters and metallicity) are known 
with errors that should be propagated to obtain the un- 
certainties on the predicted microlensing quantities (op- 
tical depth, average threshold impact parameter, number 
of events) . We have not considered these errors since they 
are far lower than the ones due to the MF uncertainties 
and the halo mass density profile. However, they should 
be correctly estimated before comparing the theoretical 
predictions with the (eventually obtained) observational 
results. On the other hand, having neglected these uncer- 
tainties does not affect our main results. 

It is interesting to compare our results to those found 
by Rhoads & Malhotra (1998) who concluded that a pixel 
lensing survey should be able to observe a total number 
of 40- 120 events after 2-5 years of monitoring of the en- 
tire globular cluster system. The method used by these 
authors is quite different from ours since they estimate 
the number of events by multiplying the event rate by 
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the effective number of fully resolved stars, defined as the 
integral of the cluster luminosity function weighted by 
a threshold impact parameter. This number is explicitly 
evaluated for NGC 7078 and then scaled to the other clus- 
ters according to the integrated absolute magnitude and 
the distance modulus. Moreover, they assume that all the 
events have the same duration (~ 30 days) , which is a very 
rough approximation, and that the halo is fully made out 
of MACHOs. If we scale our results to / = 1 and to the 
entire globular cluster system (i.e. we observe ~ 150 clus- 
ters instead of 18), we roughly get 11 - 130 events after 
three years of monitoring, which is in qualitatively good 
agreement with the results by Rhoads & Malhotra (1998), 
especially if we consider that they include also the con- 
tribution of the bulge and the disk and use different halo 
models. This simple scaling argument may be considered 
as an evidence that our results are not significatively af- 
fected by systematics induced by the method we have used 
to get them. Actually, the smaller number of events with 
respect to that predicted by Rhoads & Malhotra (1998) 
is a consequence of our most realistic assumptions on the 
halo modeling and of the lower number of targets. Thus we 
conclude that pixel lensing observations towards galactic 
globular clusters are unable to discriminate among differ- 
ent halo models because of the too small number of de- 
tectable events. 

While not useful to constrain the dark halo properties, 
globular clusters are still viable targets for a pixel lensing 
survey. In this paper, we have considered targets which 
are far above the galactic plane since we were mainly in- 
terested to the halo. However, the line of sight to many 
globular clusters goes through the disk and thus allows 
one to investigate its vertical structure and to extend the 
determination of its MF to the faint end, unreachable with 
optical observations. Furthermore, globular clusters near 
the galactic centre (such as, e.g., 47Tuc) may be used 
both as targets and to enhance the probability for lens- 
ing IjJetzer et al. 1998'jl . Then, it should be interesting to 
extend the same computations of this paper to clusters 
located towards the centre of the Milky Way or towards 
spiral arms to see whether the number of events increases. 

The analysis performed in this paper has shown that 
globular clusters are not useful targets to investigate 
the dark halo structure. However, other possible tar- 
gets are the dwarf satellites galaxies of the Milky Way 
l|Grenacher et al. 1999|l . For many of them, the line of 
sight intercepts only the dark halo and hence the (eventu- 
ally) observed microlensing events should be due to galac- 
tic MACHOs. Being dwarf galaxies dark matter domi- 
nated, it is expected that self lensing could seriously con- 
tribute to the number of events complicating the interpre- 
tation of the results. A detailed modeling of the target is 
thus needed in order to decouple the self lensing contribu- 
tion from the Miky Way dark halo one. Even if difficult, 
this task is nowaday possible (see, e.g., |Kleyna et al. 2002| 
IWilkinson et al. 20021 for Draco). In a forthcoming paper, 
we plan to repeat our analysis to see whether pixel lensing 



observations towards some dwarf satellite galaxies could 
shed light on the galactic dark halo. 
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